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^^^j^j^^GTES FOR MyLTjWAVELHg^ 
OPTICAL IMAGING 

CRQSS-REPg^gg T ° RELATEP ^"-^^ 

[0001] This application is a Continuation-In-Part 
Application of PCT/CA02/01081 , filed July 16, 2002, 
designating the United States, now pending. Thxs 
application is also a Continuation- In- Part of U.S. Patent 
Application S.N. 09/985,436, filed November 2, 2001, now 
pending. This application claims priority under 35 U.S. C. 
119(e) of U.S. Provisional Patent Application S.N. 
60/305,092 filed July 16, 2001. This application is also a 
Continuation in part of PCT/CA03/00046 , filed January 22, 
2003, designating the United States, now pending. The 
specifications of said applications are hereby incorporated 
by reference . 

MICROFICHE APPENDIX 
[0002] Not applicable. 

FIELD OF THE INVENTION 

[0003] The present invention relates to the field of 
optical imaging in which objects which diffuse light, such 
as some human body tissues, are imaged using signals 
resulting from the injection of light into the object and 
detection of the diffusion of the light in the object at a 
number of positions. More particularly, the present 
invention relates to the choice of wavelengths for 
multiwavelength optical imaging in order to provide 
enhanced information. 
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,00.4, optical -lei ae tech n iq uee 

Oo^m ana Sinuous Wave show ^ other 

for imaging breast tissue, as well the 

body parts. In TD. the objective rs to analyse t 

ooint spread function (TPSF, of an injected pulse of light 

point epreem attenuation 

as it is diffused in the tissue. With CW, 

, liaht source is measured. The information 

of a continuous light source is 

■ * tpsf and the attenuation signal is usea 

extracted from the TPSt ana 

in constructing medically useful images. 

[0005] For example, one can extract time-gated attenuation 
information from the TPSF which provides high quality 
images albeit of lower resolution than other modalities 
such as X-ray imaging. Thus, it is unclear whether the 
spatial resolution provided by optical imaging is adequate 
for diagnosing breast cancer based on morphology. 

[0006] CW and TPSF data, when processed adequately, can be 
used to extract absorption values from raw measurements. 
For example, the TPSF can be used to decouple the light 
attenuation into absorption and scattering components. This 
extra information may be clinically useful. Moreover, one 
can obtain the tissue absorption spectrum by performing 
time-domain measurements at multiple wavelengths. In tissue 
there are several molecules which absorb the light and are 
known as chromophores . Spectroscopic analysis of the tissue 
absorption spectrum permits chromophore concentrations to 
be measured. Furthermore, combination of the chromophore 
concentrations can yield physiological information, as 
opposed to morphologic information, which could provide a 
more medically useful image. 

[00071 The problem is one of knowing which are the 
dominant chromophores to include in a tissue model and then 
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-i^nt-h* to deduce their 
♦-v,^ "best" wavelengths to 
choosing the oest 

concentrations most accurately. 
SUMMARY ™ THE INVENTION 

^^7— iTc^f the invention to 
Llity * TPSP or CW-based optical images by 
SciLt combination of wavelengths and coming 
information from the combination of wavelengths. 
[0 0091 It is an object of the present invention to provide 
an objective method for choosing the wavelengths for a 
ro ultiwavelength TPSF or CW-based optical imaging approach. 
For a given set of chromophores, the best selection of the 
wavelengths is performed for the set as a whole as opposed 
to choosing the best wavelength for each chromophore 
individually. Furthermore, hardware constraints can be 
taken into consideration in order to optimize the selection 
of wavelengths for a given device. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0010] Further features and advantages of the present 
invention will become apparent from the following detailed 
description, taken in combination with the appended 
drawings, in which: 

[0011] Fig. 1 illustrates the absorption spectra used of 
oxy-Hb, deoxy-Hb, pure water and lipid; 

[0012] Fig. 2 illustrates the inverse of the condition 
number of the hemoglobin specific absorption matrix as a 
function of wavelength Xl being plotted for a system of a) 
two wavelengths, where the other wavelength is fixed at X2 = 
850 nm, b) three wavelengths where the other wavelengths 
are fixed X2 = 850 nm and X3 = 753 nm and c) four 
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waveiengths where the other wavelengths are fixed » - MO 
nm , X3 - 758 n» and M = SOO nm; 

,00,3, Pig. 3 illustrates the inverse of condition numher 
71 the specific ahsorption spectra of a*-*, and deorv- 
as a function of U and u. The plot is sy^etric w lt h 
respect to the diagonal. Regions of high values Create 
cognations of wavelengths advantageous for spectroscopy, 

,,„„, Pig. 4 illustrates the deviation of calculated 
saturation and true saturation (S (calc) -S (true) ) for a 

-1 c mm rmvn s(true) = 25%, 50% and 
model tissue containing 15 uM IHbTj , bltrue, 

75% and a lipid concentration of 40%. Two wavelengths at 
760 and B50 nm were used to fit [oxy-Hb] and [deoxy-Hb] . 
The sensitivity with respect to wrong assumptions of lipid 
and water concentrations are shown; 

[0015] Fig- 5A illustrates the inverse of condition number 
C for the specific absorption spectra of oxy-Hb and deoxy- 
Hb and lipid for a fixed wavelength X = 83 0 nm as a 
function of XI and X2 . The islands of high values indicate 
advantageous wavelengths (scaling 0-0.01); 

[0016] Fig. 5B illustrates the inverse of condition number 
C for the specific absorption spectra of oxy-Hb and deoxy- 
Hb and lipid for a fixed wavelength X3 = 830 nm as a 
function of U and X2 (same as Fig. 5A but scaling 0- 
0 . 0005) ; 

[00171 Fig. 6A illustrates the inverse of condition number 
C for the specific absorption spectra of oxy-Hb and deoxy- 
Hb, lipid and water for two fixed wavelength at X3 = 760 nm 
and X4 = 830 nm as function of Xl and X2 (scaling 0- 
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• ns of high values are advantageous for 
0.0015). Regions of mg« 

spectroscopy; 

-v,- inverse of condition number 
[001.1 Fig. 6B illustrates the averse 
c £ or the specific ahsorption spectra of *y 
Hb, lipid and water for two fixed wavelength at 
ana U - 0*0 n™ as function of U and U. ( sa.e as 
but scaling 0-0.0005); 

n -nitrates the estimation of deviations 
[0019] Fig. 1 illustrates tne 

£ ro» true saturation values for a model tissue of [HbT] 
20 S= 7S*. a lipid concentration of 40* and true wa 

concentration corresponding to 0-100, water. X «. 
wavelengths at 760. 7S0 and 8 5 0 n m were used for bac* 
calculation of S, shown here as a function of assumed water 
concentration; 

[0020J Fig. 8 illustrates the estimation of the influence 
of errors (noise) in u a on the calculated Hb concentrations 
and saturation values. A model u a -spectrum based on 20 uM 
[HbT] , £=50%, and a lipid and water concentration of 30% 
and 40% was assumed. Matrix inversion was performed for 
wavelengths 760, 730, 830 and 850 ». plotted is the change 
in calculated [oxy-Hb] , [deoxy-Hb] and saturation value 
when the u a value at a single wavelength was changed by 
,0.0001mm- 1 . This plot suggests that noise at 830 nm 
translates in the highest noise in saturation values ; 

[0021] Fig. 9 illustrates the estimation of the recovery 
of saturation values based on different wavelength 
combinations. A model tissue of 20 uM [HbT], a true 
saturation of S=75%, lipid and water concentration of 40% 
were used, in the lower plot an offset of 0.0005 mm 
independent of wavelength was added to the model tissue ua- 
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, „ offset m the upper plot) . netted are 
speetrum (no °» 6 ^ due co matrix inversion 

nations o£ the satu 

™* ^ ^ 1 , 7,0 nm and 350 nm, 2) 7,0, S30 
combinations were usea. 750 . 65 onm, 5) 

an o 050 nm. 3) 760, 780, 830 and 850 nm, 
120 -350 nm, 6) 720-900nm; and 

«m » illustrates the estimation of the recovery 

V „ values based on different wavelength 
of saturation values 

4 i i- issue of 20 MM [HbT] , a true 
combinations. A model tissue or 

* c *o% lipid and water concentration of 4CK 
saturation of S=50„ lipid ^ 

were used, in the lower plot an offset 
independent of wavelength was added to the m o del tissue ua- 
sp ectrum (no offset in the upper plot). Plotted are 
deviations of the saturation values due to matrix inversion 
and the true 75% value. The following wavelength 
combinations were used: 1) 760 nm and 850 nm, 2) 760, 630 
and 850 nm, 3) 760, 780, 830 and 850 nm, 4) 750-850nm, 5) 
720-850 nm, 6) 720-900nm. 

[00231 It will be noted that throughout the appended 
drawings, like features are identified by like reference 
numerals . 

nv.T&TIjBD DESCRIPTION OF THE INVENTION 

[0024] in accordance with the present invention, there is 
provided a method for selecting wavelengths for 
multiwavelength optical imaging. 



Tissue Chromophores 

[0025] The dominant near infrared chromophores contained 
in breast tissue are considered to be hemoglobin (Hb) m 
its oxygenated Coxy-Hb) and deoxygenated (deoxy-Hb) forms, 
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Fiq i shows the absorption spectra of 
water and Hp^s. Fig. 1 ^ 
oxy-Hh (at 10uM concentration), deoxy 

,, a r*r (100% concentration) , Hpi* 
concentration), P-e wa er U 

in " ™ * oxiaase but their absorption contr.but.on in 
cytochrome c oxidase, uui. 

tL breast is considered negligible compared to 
aforementioned chromophores . 

P hysiologic - 1 j "formation 

m Potentially useful physiological infor.atxon about 
the breast tissue can be obtained from concentrations. [] . 
of the cbromophores. The total hemoglobin concentration 
OKI. defined as tHbT] - toxy-Hb) ♦ [deoxy-Hbl , is related 

^^r,oivv cince cancer is • commonly 
to the local vascular density, bmce 

associated with an increase in vascularisation 
(angiogenesis) , a measurement of [HW] could be medically 
useful The fraction of hemoglobin that binds to oxygen is 
Known as the oxygen saturation, S. and defined as S=toxy- 
Hb] / [HbT] - increased metabolic activity increases oxygen 
demands which decreases the oxygen saturation. Since cancer 
is commonly associated with increased metabolic activity, a 
measurement of s could also be medically useful. 



Wavelength Choice 

[0027] Historically as the biomedical optics field evolved 
the wavelengths were chosen for each chromophore 
individually by observing strong near infrared spectral 
features for the given chromophore and using the closest 
hardware -available wavelength. Many researchers also used 
the isobestic wavelength of oxy-Hb and deoxy-Hb, the 
wavelength where their absorption per concentration are 
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. . chis wavelength is insensitive to the 

equal, since tins w related- to 

• .rate of the hemoglobin and can he r 
oxygenation state oi 

the CHbT] . 

the Question both posed and addressed here 
lo028] However, the J chromophores what are the 

i. that for a *>™ ^ order t0 deduce the 

optimal wavelengths to use in intere3ting to 

each chromophore? It is 
TeTt t H wavelength used by many researchers 

rrrtiot - 
- ob3ec t;r;::: ;vrr: P :L ^n, 

r:^!— . the -St action - the 

For a given wtole as opposed 

lengths is performed for the ^set chromophote 

to choosing the hest wavelength ^ 

rrrrruiie on — 3 ZTZ 

concentrations under oertain assu.pt.ons about the 
concentration,*, o £ other chromophore s, » - ~ • 
Furthermore, hardware oonstraints can also he taken , 

■ order to optimize the selection of 

consideration in order to op 

-v. f „r a aiven device. Fortunately, the recent 
wavelengths for a given a , d 

advent of turn-Key, pulsed. tunable near infrared 
length lasers has permitted more viable availability of 
near infrared wavelengths. 

,00 U One possible approach to optimize the choice of 
wavelengths for a' given set of chromophobe is to conduct ■ 
b „ te Lrce experimental study. This would cons t o 
performing numerous experiments where different 
combinations of wavelengths are evaluated for the given set 



01/16/2004 15:19 FAX 613 230 6706 



@023 



9-15186-4US-1 . 9 - 



„ s at Known eventrations until the opti^ 
of chromophobe at Known concenCra tions is found, 

cognation for deducrng be ■ htghly time - 

Obviously. this approach re UK V & ^ rf 

if is not always trivial to * 
consuming and it „ ntration9 , particularly in the 

chromophores at known concentrations, P 

case of in vivo breast tissue. 

^^^^ the serous 

POT An alternative approach which a ^ ^ 

• «f the experimental brute torus 

experiments of tne exp fol lowing 
marrix inversion sensitivity approach. While 
matrix mvei&^^ described 
example of an embodiment of the inventron wrll U de 
„ ich reference to TPSF-based parameters. r «H 

i • i i -in t-he art that tne 
• fhat bv those skilled in tne 

appreciated that Dy further 
v. ^rM-^d to CW parameters as will be furtner 
example can be adapted to ew p* 

explained below. 

l00J21 The equation which needs to be solved tan be 
written for each wavelength as: 



i 

wh ere >a is the measured absorption coefficient, ma is the 

efficient of the different 
specific absorption coefficient 

enromophores and ci is the corresponding concentration. 
[00331 This is written is matrix form as: 
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where printing in bold indicates a matrix or vector, p. is 
a vector with a number of rows corresponding to the number 
of wavelengths <*> - c is a vector with the number of rows 
corresponding to the number of chromophobe (n c ) M is a 
rectangular matrix of size n x x n c . 

[0034] If n„ = n c the system can be solved by matrix 
inversion c=M" l n a and if n x > n= the system ^is 
overdetermined and can be solved by the pseudo- inverse M* = 
(M T M)" 1 M 1 where M T is the transposed matrix of M. 

C = (M*)Ua 

[0035] The pseudoinverse M + is an n x x nc array which is 
unique. If M is square (i.e. not overdetermined) the M + = 
M' 1 . For given (i.e. chosen) wavelengths the pseudoinverse 
M + can be precalculated once and the matrix inversion 
corresponds to a simple matrix multiplication. This is the 
basis for the calculation of chromophore concentration. 

[0036] One means to quantify the expected sensitivity of a 
matrix inversion of a matrix M with respect to small errors 
in the data is the condition number C which is defined as: 

C = norm (M) • norm (M" 1 ) 

[0037] C gives an indication of the accuracy of the 
results and is an estimate of the cross-talk between the 
different channels (i.e. chromophoree concentrations). 
Values of C near 1 indicate a well-conditioned matrix, 
large values indicate an ill-conditioned matrix. The 
condition number is closely related to singular value 
decomposition (SVD) as it is the ratio of the largest and 
the smallest singular value of a matrix. 



01/16/2004 15:19 FAX 613 230 6706 



OGILVY RENAULT 



12)025 



9-15186-4US-1 _ X1 _ 

M *~>- ow-Hb and deoxy-HB at X = 760 and 
[0038] The matrix M for oxy-HD ana ^<= i 



770 nm is 



0.3871 0.1465A = 760nm 
M ^ 0.3280 0.1625A = 770nm 

[00391 A matrix inversion. is possible as the rank (N) = 2, 
however the absorption at the two wavelengths is 'similar.. 
The condition number is C = 20.49. Choosing the wavelengths 
to be X = 7 60 and 850 nm gives the matrix 

0.3871 0.1465 

M = 

0.1729 0.2645 

[0040] inspection by eye already shows that the absorption 
is very 'different'. This is confirmed by the condition 
number: C=3.206. In what follows below the inverse of the 
condition number is plotted and analyzed. It has value 
between 0 and 1. 1/C close to 1 means 'orthogonal' spectra 
and low sensitivity to cross-talk. Small values of 1/C mean 
an ill-conditioned matrix. To find the best wavelengths, 
1/C is calculated as a function of a wavelength. The 
wavelengths that give the highest values of l/C are the 
best for a calculation of chromophore concentrations and 
the subsequent physiological information such as oxygen 
saturation, S. 

[0041] Model absorption spectra were generated with the 
absorption spectra of Fig. 1 based on estimations of [HbT] , 
S, lipid and water concentration. Matrix inversion based on 
different sets of wavelengths were performed to recover 
these parameters. These parameters were compared with the 
true ones for the different wavelengths and the sensitivity 
to noise or measurement offsets considered. 
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[0042] Assuming that we fit for the hemoglobin 
concentrations only and assuming certain values for water 
and lipid concentration, for a x-wavelengths matrix 
inversion, the best combination of wavelengths to give a 
well -conditioned matrix, the sensitivity of calculated 
values of oxy-Hb and deoxy-Hb concentration and oxygen 
saturation for variations of lipid or water concentration 
and sensitivity of S to measurement noise have been 
determined. 

[0043] In Fig. 2 the inverse of the condition number is 
shown for matrices of oxy-Hb and deoxy-Hb specific 
absorption coefficients for 2, 3 and 4 wavelengths. In each 
case one wavelength (X x ) was varied between 650 and 950 nm 
while the remaining wavelengths were fixed X 2 = B50 nm (2- 
wavelength system) , X 2 = 650 nm and X 3 = 758 nm (3- 
wavelength system) , and X? = 850 nm, X 3 = 758 nm and X 4 = 
800 nm (4-wavelength system). Fig. 2 indicates that the 
selection of two wavelength at X 1 = 850 nm and k 2 = 700 nm 
gives the highest values of 1/C and when the wavelength 
range is restricted via hardware constraints to > 750 nm, a 
system that includes the peak wavelength of deoxy-Hb close 
to 760 nm is advantageous. It does not matter whether two 
or more wavelengths are used. This somewhat 
counterintuitive result is valid only without measurement 
noise and noise in the background absorption. 

[0044] Fig. 3 further highlights this finding for a two- 
wavelength matrix inversion. In this figure l/C is plotted 
as a function of both at Xi and X 2 in the range 650-950 nm. 
The plot is symmetric with respect to the diagonal. Regions 
of high 1/C -values can be chosen and the corresponding 
•good* wavelengths can be read off the axis. It is apparent 
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that (with the restriction to > 750 nm) the one wavelength 
should be close to 760 nm while the other one can be within 
the range 830-900 nm without substantially affecting the 
condition number. 

[0045] Using the spectra shown in Fig.l, model tissue 
absorption spectra were generated. Based on matrix 
inversion values of [oxy-Hb] , [deoxy-Hb] and S were 
backcalculated and the sensitivity to incorrect assumptions 
about the [water] and [lipid] tested. One approach is to 
take the measured u a spectra and subtract water and lipid 
absorption corresponding to certain assumed concentrations. 
For the data shown in Fig. 4, a model tissue containing 15 
\M [HbT] , (true) saturation values of S = 25%, 50% and 75% 
was used. Lipid concentration was 40%. It was tested how a 
misjudgement of water concentration affects the 
recalculated S value. To test the error in a simple two- 
wavelengths-fit (760 and 850 nm) , the assumed lipid 
concentration was varied between 0 and 100%. When the 
assumed water concentration is right (lower three lines in 
Fig. 4) , the deviation in saturation between true and 
calculated values is < +2% (obviously with zero error for 
the right lipid concentration of 40%) . A misjudgement about 
the water concentration by 20% (upper lines in Fig. 5) 
results in additional errors in S (calc) -S (true) of up to 2% 
for S=75%, 4% for S=50% and 8% for S=25%. These errors in S 
are a function of the underlying tissue absorption 
coefficients. The values here give an indication about the 
order of magnitude. 

[0046] Having a system with more than two fit-parameters, 
best wavelength combinations, for a three- components system 
of oxy-Hb, deoxy-Hb and lipid system, for a four- components 
system of oxy-Hb, deoxy-Hb, lipid and water, and the 
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sensitivity of calcuiation of s to noise at the different 
wavelengths have been determined. 

[0047] in Figs. 5A and SB the inverse of the condition 
number is plotted for a three wavelengths system based on 
the oxy-Hb, deoxy-Hb and lipid specific absorption spectra 
as a function of X, and X 2 . The third wavelength was fixed 
at X3 = 830 nm. Again, the plot is symmetric with respect 
to the diagonal. From Fig. 5A it is apparent, that there 
are three "islands" of high 1/C values. Unfortunately, all 
of these island would include wavelengths outside an 
imposed hardware constrained wavelength range of 750 to 850 
nm. Plotting the same data in a different scale (Fig. 5B) 
shows that there is just a single preferential combination 
within this hardware constrained wavelength range: 760 and 
780 nm. 

[0048] Equivalent to Figs. 5A and 5B, the inverse of C for 
a 4-wavelengths system is plotted in Fig. 6A and 6B. Again, 
the difference between them is the scaling. Two wavelengths 
were fixed at X3 = 760 nm and X4 = 830 nm. Including the 
wavelengths outside the 750-850 nm range there appear four 
preferential combinations. Restricting the wavelength range 
to 750-850 nm there are just two advantageous region 
(marked by the white rectangle in Fig. 6B) : 780 nm and 850 
nm, and 780 nm and 815 nm. 

[0049] From the analysis based on matrix condition 
numbers, the best wavelength combinations for 2, 3 and 4 
wavelengths measurements are the following: 
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Table 1 



Wavelength combinations for 2, 3 and 4 wavelengths measurements 





Wavelength 
Range 


Best wavelengths (nm) 
M X2 X3 X4 


fit for 


see 
Fig. 


2-X 


650-950 nm 
750-850 nm 


700 >860 
760 850 


Oxy-Hb. deoxy-Hb 


2,3 


z-x 


650-950 nm 
750-850 nm 


700-760830 925 
830 860-870 925 
760 780 830 


+ lipid 


A,B 


4-X 


650-950 nm 
750-850 nm 


760 830 860 925 
700 760 830 925 
760 780 830 850 
760 780 815 830 


+ lipid, + water 
Best combination 


6A, 
B 



[0050] Furthermore, it must be pointed out that including 
more wavelengths doea not increase the condition number. 
E.g. for the four chromophores and all wavelengths in the 
range 750-850 nm, l/C = 0.000314. This is lower than the 
value (l/C = 0.00036, compare with Fig. 6B) when only four 
wavelengths (.760, 7 60, 83 0 and 650 nm) are used. In a 
system without noise and no other chromophores than the 
four considered here, a 4 -wavelengths system is the 
optimal . 

[0051] While certainly only a 4 -wavelengths measurement 
allows [oxy-Hb] , [deoxy-Hb] , [lipid] and [water] to be 
determined, and a 2 -wavelengths system (see Fig. 5) is not 
sufficient, the question is posed whether a 3 -wavelengths 
measurement might supply S values with a high enough 
precision. In this case the concentration of one 
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chromophore (water or lipid) must be guessed and the 
corresponding absorption subtracted from the measured u a 
values. This was tested with a model absorption spectrum 
and is shown in Pig. 7 for 760, 780 and 850 nm. True water 
concentration was varied between 10 and 100% (the different 
lines), and the difference between calculated and true 
saturation values plotted as a function of assumed water 
concentration. For instance, for a true water concentration 
of 50%, a misjudgment of the water concentration by 10% 
results in an error in S by about 4%. 

[0052] Up to now, only "perfect" data sets were considered 
with no noise. In real situations there are problems due to 
measurement noise that is random for the different 
wavelengths; unknown chromophores in the tissue, i.e. there 
is a background absorption coefficient the spectrum of 
which we do not know; and possible systematic errors in the 
primary u a recovery. 

[0053] There are an ample number of parameters which can 
be considered and as examples the following two questions 
are considered. First, is the oxygen saturation more 
susceptible to noise at certain wavelengths? Second, what 
is the influence of an offset in the u a data? 

[0054] The influence of errors (noise) in u 6 on the 
calculated Hb concentrations and saturation values was 
estimated in a model tissue based on 20m [HbT] , S= 50% and 
a lipid and water concentration of 30% and 40% 
respectively. Matrix inversion was performed on the u a - 
spectrum of this model tissue for wavelengths 760, 780, 830 
and 850 nm. The sensitivity to noise (i.e. variations in 
jia) at the different wavelengths wa3 estimated by changing 
the absorption coefficient at a single wavelength by 
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the change in calculated loxy-Hb] , 
+0.0001 mm' 1 , in Fig. 8. the cnang ,„ Ql ... 
taeoxy-Hb] and oxygen saturation value due to this noise 
plotted. This «gure shows that the change in oxygen 
saturation value is aoout -2% for changes at 7,0 nm «.£ 
at 700nm. while it translates to a variation of at 330 



nm. 



l0 055] There are two conditions that produce an offset m 
the measured u 8 -spectra with respect to the true values. 
First, the algorithm for ^-calculation based on TPSF-based 
optical imaging might lead to a systematic offset e.g. due 
to residual crosstalk between absorption and scattering 
parameters. Second, the tissue absorption might have a 
background of unknown origin (chromophore) . Under both 
conditions the fitting of u a -data with the four 
chromophores is hampered. The effect of such an offset for 
different wavelength combinations is estimated with a model 
spectrum of 20 pM [HbT] , S=75% and water and lipid 
concentration of 40%. An offset of 0.0005 mm" 1 was added to 
the p.- values independent of wavelength. The effect on the 
calculated oxygen saturation values is known in Fig. 9 for 
combination of 2, 3 and 4 wavelengths as well as continuous 
spectra between 750-850, 720-850 and 720-900 nm. It is 
apparent that the lowest error in S is achieved by the 4- 
wavelengths combination. Including more wavelengths 
increases the error. In Fig. 10 the same calculation was 
done, however, for a true oxygen saturation value of 50%. 
Here the lowest error is achieved by the 720-850 nm 
wavelength range, while using less wavelengths or 
increasing the fitting range to 900 nm results in larger 
errors . 
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[0056] Based on che assumption that the dominant tissue 
chromophores are oxy-Hb, deoxy-Hb, water and lipid and 
analysis of the matrix condition number, measurements at 
the wavelengths 760, 780, 830 and 350 nm supply an optimal 
data set when the wavelength range is limited to 750-350 nm 
under ideal conditions. As shown in Figs. 5 and 6, 
inclusion of shorter and longer wavelengths promise a 
better matrix inversion. Under real conditions there is .no 
clear-cut answer about the improvement when more 
wavelengths are included (see Figs. 9 and 10). It might be 
advantageous to reduce measurement noise at 4 wavelengths 
due to longer scan times rather than to include more 
wavelengths. As demonstrated in Fig. 8, to achieve an 
optimal accuracy the noise level at different wavelengths 
has to be adjusted which might require different 
measurement times at certain wavelengths . 

[0057] Strictly speaking the work presented here was 
achieved by optimizing a 2-wavelength system and then 
optimizing a 4 -wavelength system where 2 of the wavelengths 
were fixed at the optimized 2-wavelength system values. 
Whilst this is easier to display graphically, preferably 
all 4 wavelengths would be permitted to vary in a global 
optimization process. Fortunately, for the specific example 
presented here when all 4 wavelengths are permitted to vary 
the same optimal solution is found. However, this may not 
be true for all situations and a global optimization is 
preferred. 

[0058] It will be appreciated that parameters of the 
system other than wavelengths can be optimized. For example 
the light source type and source power, the detector type 
and detector aperture for each wavelengths, the choice of 
image algorithm, the source/detector geometries, the 
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acquisition time and the noise characteristics are 
parameters that can be adjusted or chosen, as would be 
known to one skilled in the art, to optimize the 
determination of chromophobe concentration and the optical 
image obtained therewith. In a preferred embodiment, 
optimization of the parameters is performed by taking in 
consideration the matrix inversion approach for optimizing 
wavelength selection as described above. 

[0059] It is also understood that it will be obvious to 
those skilled in the art that the same approach for 
choosing optimal wavelengths can be applied to optical 
absorption spectroscopy in general. For example, in other 
embodiments of the present invention the method of the 
present invention is also used for choosing the optimal 
wavelengths for analyzing the components of paints, 
pharmaceutical products, food, grain or any other turbid 
media. 

[0060] It is also understood that the proposed method 
applies both to the analysis of absolute chromophore 
concentrations as to their changes or relative 
concentrations . 

[0061] It is also understood that the proposed method 
applies both to the analysis of absolute chromophore 
concentrations as to their changes or relative 
concentrations . 

[0062] Further, it is also understood that the proposed 
method applies when continuous wave (CW) methods are used 
and an assumption for scattering is made, e.g. a constant 
value or following a scatter-power wavelength dependent law 
, in order to infer the absorption coefficient, racher than 



01/16/2004 15:22 FAX 613 230 6706 



OGILVY RENAULT 



@034 



9-15186-4US-1 _ 2Q _ 

measuring the absorption coefficient directly with a TPSF- 
based approach. 

[0063] Thus, changes in CW measurements can be converted 
into absorption changes M(A) which are then used to 
calculate changes in concentrations*, based on a modified 
Beer-Lambert law which assumes a predetermined wavelength 
dependence of the optical pathlength D.W (Cope * Delpy, 
1998; Bssenpreis et al.. 1932) to account for scattering. 

[0064] Thus Ac i =N" 1 AA / where N is similar to M but 
incorporates DM) • and is the matrix which is used for 
optimizing the choice of wavelengths. 

[0065] Further, it is understood that the proposed method 
applies not only for a mixture of endogenous absorbers 
(chromophores) , but also for mixtures of exogenous 
absorbers with known spectra, such as optical dyes or 
fluorophores, or for a mixture of both endogenous and 
exogenous absorbers. 

[0066] While the invention has been described in 
connection with specific embodiments thereof, it will be 
understood that it is capable of further modifications and 
this application is intended to cover any variations, uses, 
or adaptations of the invention following, in general, the 
principles of the invention and including such departures 
from the present disclosures as come within known or 
customary practice within the art to which the invention 
pertains and as may be applied to the essential features 
herein before set forth, and as follows in the scope of the 
appended claims. 
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l00 671 The embodiment (s) of the invention ascribed above 
is (are) intended to be exemplary only. The scope of the 
invention is therefore intended to be limited solely by the 
scope of the appended claims. 



